% Muenchow flower (waiting time) data 
% female is a 1x47 vector of waiting times
% cen is a 1x47 vector or censor indicators
% flowers is a 47x2 vector
% s0 = kaplan meier estimate
female=[1 2 4 4 5 6 7 7 8 8 8 9 14 15 18 18 19 23 23 26 28 29 29 29 30 32 35 35 37 39 43 56 57 59 67 71 75 75 78 81 90 94 96 96 100 102 105];
censor=[0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 1 0 1 1 0 1 1 1 1];
cen=1-censor;
flowers=[female',cen'];
s0=KMcdfSM(female',cen',1);